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Abstract 



The extraction efficiency of ultracold neutrons from cryogenic 
moderators depends critically on the neutron transparency 
of the moderator material. A Monte Carlo simulation of the 
probability of the UCN going out from non-ideal (disordered) 
sohd deuterium crystals has been performed. It was based on 
the use of a correlation function describing the density fluctu- 
ations in a disordered material, the latter being inferred from 
the measured very low neutron energy total cross sections for 
this material. 



Ultracold neutrons [T] (UCN, the energy below ~ 0.3/ieV), which can be con- 
tained in closed volumes, have found important applications in the investigation 
of fundamental properties of the neutron [2]. 

The best UCN source - ILL UCN turbine [3] yields maximum UCN density in 
a trap of ~ 50 n/cm^ and, respectively, the flux ~ 5 x lO^n cm'^s^^. 

Recent years brought hope for significant progress in the intensity of very low 
energy neutron beams. This is connected with the possible use of the most effective 
cold moderators: solid deuterium or solid deuterocarbons at low temperatures of 
~ 5 - 10 K for production of very cold neutrons. 

The most promising and popular material for the UCN converter is solid deu- 
terium. The first test of the UCN production in solid deuterium has been per- 
formed by the PNPI group jl] , which then has been studying experimentally the 
UCN production in deuterium during many years [5-11]. 

The first estimate of the UCN production and up-scattering cross sections 
in solid deuterium was made in [12]. The phonon frequency spectrum of solid 
deuterium was calculated in [13] using the dispersion relations measured in [TH . 



Later calculations in the incoherent approximation [TH] of the UCN scattering on 
phonons in solid deuterium reproduced the results of [13], this publication also 
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contains the calculation of the UCN spin-flip up-scattering on para-deuterium with 
para-ortho transitions of deuterium molecules. The incoherent approximation and 
the density of states of [I^ was used in the calculations of the UCN production 
in |l6j. 

It was proposed that because of comparatively short neutron lifetime in solid 
deuterium (140 ms in pure ortho-deuterium at T=0 K) solid deuterium can find 
better application in a pulsed mode of the UCN production [17,9]. 

A solid deuterium UCN source was also demonstrated in the experiments in 
LANL [TS], it is under construction at FRM-2 jTD] and is in the stage of completion 
in PSI [20]. 

The incoherent model of neutron interaction in solid deuterium was generally 
confirmed in the later works ITHl EH 1221 123, 



Recently a re-evaluation of the neutron scattering in solid deuterium was per- 
formed in 1251, it based on the calculations in the frame of a rigorous model 
of coherent neutron scattering by deuterium nuclei. The resulting scattering cross 
sections obtained in this work are significantly different from the previous pre- 
dictions 1121 [131 [15] based on the incoherent approximation. Especially large 
difference between the incoherent approach and coherent one was found in the 
case of UCN up-scattering with phonon annihilation in the temperature range of 
5-10 K: the first one 3-4 times surpassed the second one. These new calculations 
found some experimental confirmation in [i26]. 

It was clear from the very beginning that structural homogeneity of a deu- 
terium crystal is very important for the UCN yield from the UCN converter. 
The neutron scattering on structural inhomogeneities increases the time neutrons 
spend inside the UCN converter before leaving it, increasing the probability of 
their capture and up-scattering. It was discovered in |[27j that the transmission 
of very low energy neutrons through solid deuterium depends critically on the 
way of preparation of a solid deuterium sample. The light transmission investi- 
gations of solid deuterium crystals were performed in |j2E]- It was demonstrated 
that optical quality of crystals depends on the speed of crystal growing - more 
transparent crystals were obtained when they were carefully grown over a time 
period of about 12 hours at temperature close to the triple point. 

The total cross section in the UCN energy range for solid deuterium frozen 
from liquid state was measured in [22]. These measurements showed that the 
best transmission have the solid deuterium samples prepared in the process of 
very slow cooling - less than 1 K per hour. This procedure of obtaining good qual- 
ity deuterium crystals was based on preliminary investigations described in fIE\ . 
Temperature cycling between 5 K and 10 K, and between 5 K and 18 K seriously 



deteriorated the neutron transmission. It is obvious that crystal imperfections 
introduced by thermal cycling increased the UCN scattering. In contrast to this 
experience, the solid deuterium crystals used as UCN converters and condensed 
directly from gas phase significantly increased the UCN yield being subjected to 
repeated thermal cycling \22\ . 

We consider solid deuterium crystals frozen in the process of very slow cooling 
[22] as the most perfect of all possible samples prepared from liquid, and use the 



m 



energy dependence of the total cross section (Fig. 2 of [29]) to infer the parameters 
characterizing structural inhomogeneity of the sample. These parameters are 
used further to simulate the probability for the UCN generated in a similar solid 
deuterium UCN source, to go out from the source to vacuum. 

To obtain the UCN scattering cross section on structural inhomogeneities the 
experimental points of Fig. 2 of Ref. [l29j were corrected for the incoherent elastic 
scattering on deuterium (2.05 6), the UCN capture by deuterium (1.19 b/v) and 
hydrogen (733 b/v) nuclei, the UCN up-scattering in the result of one-phonon 
annihilation (0.83 b/v), and the spin-flip up-scattering on para-deuterium (123 
b/v), where v is the neutron velocity in medium in m/s. The cross sections of 
the latter processes were taken from |[T5]. All the microscopic cross sections here 
and further are given per one atom of respective element. The concentrations 
of hydrogen (0.05%) and para-fraction (1.4%) were taken as they were done in 
p9j . The atomic density of solid deuterium at 5 K was taken to be 6 x 10^^ 
cm~^. The resulting neutron wavelength dependence of the UCN scattering on 
inhomogeneities is shown in Fig. 1. 

The transmission of very slow neutrons through an inhomogeneous medium was 
studied previously by A. Steyerl [30]. It was demonstrated that transmission as a 
function of neutron wave length may be used to deduce characteristic parameters 
of inhomogeneities, in particular their size and density. 

It is known that in the Born approximation the differential macroscopic cross 
section of elastic scattering for neutrons transmitting through an isotropic inho- 
mogeneous medium has the form 



- G{p) p dp, (1^ 



where m is the neutron mass, q=|/c' — /c| is the neutron wave vector change, and 
G(r, r')=< (5U(r) ^U(?^) >, {p = |?^ — ?^) is the correlation function of fluctuations 
of the local neutron-medium interaction potential. 
The latter is 



where Nj is the atomic density and bj is the coherent scattering lengths of nuclei 
of the medium, so that (5U(r)=U(r) - <U(r) >. 

The total cross section after integration over solid angle with the account of 
the solid angle of the neutron detector is: 



(3) 



CTTL \ ^ 1 roo 
^ ) p Jo G{p)[cos{2kp ■ sinOo) - cos{2kp)]dp, 

where 26q is the angle of the neutron detector from the sample. 
For the exponential correlation function 

G{p) = Goe-'IP\ (4) 

where po is the correlation length, we have the expression 



for the total cross section, and 



(6) 



for the differential cross section. 

This formalism was used for the interpretation of the measurements of the total 
cross sections for solid deuterium at 5 K [29] . The parameters of the correlation 
function of Eqs. (4) and (5) were inferred by the least squared method from the 
corrected experimental points for the macroscopic cross section of the UCN scat- 
tering on inhomogeneities of Fig. 1. The inferred parameters of the exponential 
correlation function were found to be: Go = (0.33 ± 0.08) neV^, pQ = (28 ± 3.3) A. 
The total cross section of Eq. (5) with these parameters is shown in Fig. 1. 

Fig. 2 shows the macroscopic total and transport cross sections of scattering 
on inhomogeneities: 

^tr = /(I - cosO)d^, (7) 

where 9 is the neutron scattering angle. 

The difference between total and transport cross sections is about 5% at the 
neutron energy of 200 neV and is lower at lower energies, it follows that the 
scattering by inhomogeneities is almost isotropic in this energy range. 

The differential cross section of Eq. (6) with inferred parameters Go and po 
was used to perform Monte Carlo simulation of the UCN transport in a cylindrical 
deuterium crystal aiming to obtain the probability for the UCN to go out of the 
crystal. 



The Monte Carlo program used here is the development of our program used 
earlier for simulation of transport of very slow neutrons in matter, in particular 
in cold moderators with their temperature changing during UCN motion in a 
moderator [32j . 

The UCN were assumed to be born homogeneously with isotropic angular 
distribution, the reflection from the side walls of a cylinder was specular, from 
the top and bottom walls - according to the cosine law. The quantum reflection 
from the boundary deuterium-vacuum at the exit side of the sample was taken 
into account according to: 

R=('-^^^^)\ (8) 

where R is the neutron reflection probability, vj_ is the normal to the boundary 
component of the neutron velocity in medium, vo^± = (-y^^ + 'U^)^/^ is the normal to 
the boundary component of the neutron velocity in vacuum, Vb=4:.4:6 m/s is the 
calculated boundary velocity of deuterium, the boundary energy Eb = 104 neV. 
Recent measurement [|33j of the minimal neutron energy from the solid deuterium 
source Emin = (99 ± 7) neV, which should correspond to the boundary energy of 
deuterium, does not contradict to this value. 

Two variants for the thickness of the solid deuterium sample were chosen: 15 
cm - mean thickness of the solid deuterium source in PSI [l20j and 8 cm - typical 
for the UCN source at the TRIGA reactor of Mainz University |ISl 122]. 

The UCN up-scattering due to phonon annihilation was calculated in two vari- 
ants: the incoherent model jlS] and coherent one [l25j. The result of this simula- 
tion as the probability for the UCN to go out from the cylindrical solid deuterium 
crystal are shown in Figs. (3-8). 

They illustrate effects of temperature, concentration of para-fraction in deu- 
terium, and inhomogeneity of solid deuterium on the probability for the neutrons 
to leave the moderator. 

Fig. 3 shows the probability for the UCN to go out from perfectly homogeneous 
solid deuterium cylinder with a thickness of 15 cm at different temperatures, at the 
concentration of para-fraction Cp = 0, the up-scattering on phonons was calculated 
according to the incoherent model of Ref. [TS] and the coherent model of Ref. 
[ i25j . The following one-phonon annihilation cross sections were used: for the 
incoherent variant 0.83 b/v, 2.78 b/v, and 7.9 b/v for temperatures 5, 7, and 
9 K, respectively; for the coherent variant 0.2 b/v, 0.87 b/v, and 2.78 b/v for 
the same temperatures. Further the effect of phonon upscattering on the UCN 
yield will be discussed only for the variant of calculations in the incoherent model, 
the calculations according to the coherent model predict much lower effect of the 



phonon annihilation on the UCN loss in a moderator. It is seen, that the phonon 
upscattering effect of heating of an ideal deuterium crystal from 5 to 9 K decreases 
the UCN yield more than 3 times down to ~6.5% of produced neutrons with an 
energy 150 neV. 

Fig. 4 shows the same as in Fig. 3 (thickness of 15 cm, Cp = 0), but for the 
inhomogeneous solid deuterium with the inferred parameters of the exponential 
correlation function Go = 0.33 neV^, po = 28 A. The effect of inhomogeneities 
is not very essential for this rather good crystal: neutron yield is decreased 1.3 
times and 1.2 times at temperatures 5 and 9 K, respectively, compared to an ideal 
crystal for the same UCN energy. 

The effect of spin-flip upscattering due to admixture of para-deuterium in an in- 
homogeneous crystal with a thickness of 15 cm at 9 K is seen in Fig. 5. Compared 
to the effect of phonon annihilation, additional influence of the UCN upscattering 
with para-ortho transition of deuterium molecules is not as dangerous: at 2% of 
para-deuterium the UCN yield decreases 22% for neutrons at 150 neV. 

Naturally, all the neutron loss effects are larger for thicker sample. Figures 
6, 7, and 8 show the neutron energy dependence of the effects of temperature, 
concentration of para-fraction in deuterium, and inhomogeneity of the crystal on 
the probability for the neutrons to leave the deuterium moderator with a thickness 
of 8 cm. 

Fig. 6 shows the probability for the UCN to go out from a perfectly homoge- 
neous solid deuterium cylinder with a thickness of 8 cm at different temperatures, 
at the concentration of para-fraction Cp = 0. In this case heating of an ideal 
deuterium crystal from 5 to 9 K decreases the UCN yield due to the phonon- 
upscattering effect the same 3 times, but the absolute values of the probability 
for the neutrons at 150 neV to leave the moderator are larger: 0.36 and 0.12 
compared to 0.21 and 0.065 at these temperatures at a thickness of 15 cm. 

Fig. 7 shows the same as in Fig. 4 (inhomogeneous solid deuterium, = 0), 
but for a thickness of 8 cm. The effect of inhomogeneities is not essential in this 
case too: neutron yield is decreased 1.24 times and 1.2 times at temperatures 5 
and 9 K, respectively, compared to an ideal crystal for the same UCN energy. 

And, at last. Fig. 8 shows the effect of spin-flip upscattering due to admixture 
of para-deuterium in inhomogeneous crystal with a thickness of 8 cm at 9 K. 
Compared to the phonon annihilation additional UCN upscattering with para- 
ortho transition of the deuterium molecules is not as dangerous: at a concentration 
of para-fraction 2% the UCN yield decreases approximately the same 20% for 
neutrons at 150 neV. 

Maximum UCN losses are seen if we compare, for example, the probability for 



the neutrons to leave the sample at the ideal condition of homogeneous material 
at T=0 and in the worst shown variant: inhomogeneous solid deuterium at T=9 
K, concentration of para-fraction Cp = 0.02. In this case for the UCN at 150 neV 
the yield is 6.5 times lower than from the ideal crystal with a thickness of 15 cm 
and is only 4.5% of produced neutrons. At temperature of 5 K the UCN yield 
from a homogeneous crystal is 4.2-4.7 times larger than from a inhomogeneous 
one at 9 K depending on the para-fraction concentration. For deuterium thickness 
8 cm maximum overall losses are lower than for 15 cm: the yield is about 8% of 
produced neutrons at 9 K and 2% concentration of para-deuterium molecules. 

The author is grateful to Dr. Malgorzata Kasprzak for sending the data tables 
for Fig. 2 of Ref. [l29j. 
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Figure 1: Cross section of scattering on inhomogeneities for the deuterium sample at 5 K (Fig. 
2 of [29]) as a function of the in-medium neutron wave length after subtraction of the neutron 
capture, up-scattering, and the elastic incoherent scattering by deuterium nuclei. The curve is the 
approximation according to Eq. (5) with the deduced parameters (^0=0. 33 neV^, po=28A. 




Figure 2: Macroscopic total and transport cross sections of neutron scattering by the solid deu- 
terium crystal inhomogeneities according to Eqs. (5) and (6) with the deduced parameters 
^0=0.33 neV2, po=28A. 
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Figure 3: The probability for the UCN to go out from the homogeneous sohd deuterium cyhnder 
with thickness 15 cm, the concentration of para-fraction Cp = 0, the up-scattering on phonons was 
calculated according to the coherent model of Ref. [25] and the incoherent model of Ref. [T5] . 



w^^j D-inhomog.,d=15cm, c=0 




neutron energy in vacuum [neV] 



Figure 4: The same as in Fig. 3, but for the inhomogeneous sohd deuterium with the inferred 
parameters of exponential correlation function Gq = 0.33 neV^, po = 28A 
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Figure 5: The probability for the UCN to go out from the inhomogeneous sohd deuterium crystal 
with thickness 15 cm and temperature T=9 K, at two concentrations of para-fraction Cp = 0.01 
and 0.02, the up-scattering on phonons was calculated according to the coherent model of Ref. 
[25] and the incoherent model of Ref. jl5j . 




Figure 6: The same as in Fig. 3, but for the deuterium thickness 8 cm. 




Figure 7: The same as in Fig. 4, but for the deuterium thickness 8 cm. 
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Figure 8: The same as in Fig. 5, but for the deuterium thickness 8 cm. 



